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Abstract. Focusing attention on top-quark pair production and its decay processes at tiie LHC, non-standard 
top-quark couplings are studied based on tlie effective Lagrangian constructed with 5 L'(3)x5 t/(2)xL'(l) invari- 
ant dimension-6 operators. The optimal-observable analysis is carried out for the charged-lepton distributions 
in pp — > tiX — > £*X' (f = e or jS) in order to estimate the expected statistical uncertainties in measurements of 
those non-standard ttg and tbW couplings that contribute to this process in the leading order. 
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1 Introduction 

Almost twenty years have passed since the top quark 
was discovered at the Tevatron [T, Even today, how- 
ever, this quark is still the heaviest in the experimentally - 
confirmed elementary particles. Therefore, precise mea- 
surements of the top quark are expected to be one of the 
most promising ways for finding possible signals from 
new physics beyond the standard model in the current sit- 
uation where there has been no discovery of non-standard 
particles 1 3 ] . The Large Hadron Collider (LHC) is now the 
unique facility for those top-quark studies in stead of the 
Tevatron. 

In any processes, non-standard signals mean non- 
negligible differences between actual experimental data 
and the corresponding standard-model prediction. Such 
deviations might be originated from quantum effects of 
new particles which are not belonging to the standard- 
model framework. Those quantum effects are usually 
characterized as new form factors corresponding to the 
non-standard couplings. The size of the form factors are 
calculable if you take a specific model which is an exten- 
sion of the standard model. We are also able to treat those 
form factors as parameters in the manner of effective- 
Lagrangian approach, if we would like to perform anal- 
yses as model-independently as possible. The latter is the 
strategy which we have been taking for our new-physics 
search in top-couplings |4-6]. 

In this note, we focus on our latest work |6] and show 
the studies on the top-quark production and its semilep- 
tonic decay process pp — » ttX — > {^X' ({ = e ot fi) 
at the LHC aiming to estimate the statistical significance 
of measurable non-standard top-quark couplings. Here, 
non-standard interactions which are originated from the 



quantum effects of non-standard particles are described by 
SU(3) X SU(2) X t/(l) invariant dimension-6 operators, 
and parameterized as non-standard top couplings. The 
statistical significances are estimated through the optimal- 
observable analysis. 

2 Strategy 

In this section, we explain briefly the effective Lagrangian 
and the optimal-observable analysis, both of which play 
important roles in our analysis. 

2.1 Effective Lagrangian 

The effective-Lagrangian approach, in which its low- 
energy form reproduces the standard model, is one of the 
most general methods to describe new-physics phenom- 
ena when the energy of our experimental facility is not 
high enough to produce new particles. Assuming any non- 
standard-model particles too heavy to appear as real ones, 
the effective Lagrangian is given as 



Xeff = + ^ J] ( CiOi + c;ol ) , 



(1) 



where Xsm is the standard-model Lagrangian, O, mean 
S U(3) X S U(2) X t/(l) gauge-invariant operators of mass- 
dimension 6 involving only the standard-model fields ^ 
|8|] and their coefficients C,- parameterize quantum effects 
of new particles at an energy less than the assumed new- 
physics scale aQ In this framework, all the form factors 
related to C, are dealt with as constant parameters, with- 
out supposing any specific new-physics models. 



^e-mail: hioki@ias.tokushima-u.ac.jp 
''e-maik ohkuma@fukui-ut.ac.jp 



'since dimension-5 operators violate the lepton-number conserva- 
tion, they are not treated hereafter. Therefore, we deal with dimension-6 
operators, which give the largest contributions in relevant processes. 



A full list of the dimension-6 operators were initially 
presented by Buchmiiller and Wyler J?!, and also by Arzt 
et al. fs']. However, it was pointed out that some opera- 
tors there are related with others through equations of mo- 
tion, which means that they are not independent of each 
other |9|]. Then, the whole related operators were rear- 
ranged to get rid of this redundancy in Refs. lllOlll ill . 

Following the notation of Ref. ifioll . the relevant ef- 
fective Lagrangian which describes non-standard interac- 
tions of the third generation for the parton-level process 



qqigg ■ 



is given in [5J as 



-Ceff - -Ctig.gg + -ClbW '■ (2) 



- iArWr {dv + idAy5Wx)G%(x) 1, (3) 

+ Mx)—(f2PL + f^PR)Mx)d^w;(x) ], 



(4) 



where gs and g are the S U (3) and S U (2) coupling con- 
stants. Pi/ It = (1 + 75)/2, dv, dfy and f^f are form factors 
defined as 



_ V2WM, 33 _ V2WM, 33 



2 2 

= _ ^/2r^^*— f « = - V2r^^ — 



(5) 



with u being the Higgs vacuum expectation value and Vtb 
being (tb) element of Kobayashi-Maskawa matrix. In par- 
ticular, dv and c/^ are the so-called chromomagnetic- and 
chromoelectric-dipole moments, respectively. 

In the following calculations, we use the above ef- 
fective Lagrangian for top-quark interactions and the 
usual standard-model Lagrangian for the other interactions 
which are not affected by top quarks hereafter. 

2.2 Optimal-observable analysis 

The optimal-observable analysis is a method that could 
systematically estimate the expected statistical uncertain- 
ties of measurable parameters lll2l - [l5ll . Here we give a 
brief review of this procedure. 

First, assume that we have a cross section 



dcr 



(6) 



where are known functions of the final-state variables 
and c, 's are model-dependent coefficients. The goal is to 
determine the c,'s. This can be done by using appropriate 
weighting functions Wi{^) such that J Wi((p)2J((p)d(f> - c-,. 
Then, we determine as the next step. In general, dif- 
ferent choices for are possible, but there is a unique 



choice for which the resultant statistical error is mini- 
mized. Such functions are given by 



(7) 



where Xij is the inverse matrix of Mij which is defined as 
Mij = j d<f>M4,)fj(cf>)/I(<p). (8) 

Finally, using Eqs. O and ([8]), the statistical uncertainty 
of c, is obtained as 



\Sci\ = ^JX,(rT/N, 



(9) 



where a-f = J (dcr / d(/))d(/) and is the total number of 
events. 

Let us apply this technique to the process pp ttX 
{^X', expressing its differential cross section (the angular 
and energy distribution of the charged-lepton £'^) which 
we derived in [5J as follows: 



d^CTf 



dEcd cos 6c 



= fsuiEc, cos 0t) + dvfdyiEf , cos 0() 

+ dRfd„(E(, cos Oc) + dlf^iiEt, cos 0e) 
+ dlfAEe, cos Ot) + ■■■ , (10) 



where fsuiEf, cos 6[) denotes the standard-model con- 
tribution, all the other fi{Ec,cos6c) describe the non- 
standard model terms corresponding to their coefficients, 
and du is defined as 



dR . ^Re(/«). 

m, 



(11) 



The explicit forms of fi{Ec, cos Be) at the parton level are 
easily found in the relevant formulas in 

Since the magnitude of dv and dA has been shown 
small Jsl |6t], we neglect any contribution from terms 
quadratic (or higher) in those non-standard model param- 
eters hereaftero Note that all dA contributions disappear 
under this linear approximation, because there is no term 
proportional to dA- Thus, the angular and energy distribu- 
tions of a decay lepton are written as 



d cos Of 



- gi(cos0() + dv giicosOe), 



(12) 



5 = hi(Ee) + dv h2(E() + dR h3(E[), (13) 
dE[ 

where giicos 6c) and hi(Ec) are given by 

ef/(cos 6c) = J dEc fi(Ee, cos 6c), 



hi(Ec) 



I 



d cos 6c fi(Ec, cos 6c) 



with / - 1,2 and 3 corresponding to / = SM, dy and 
dR, respectively. Here should be one comment about the 



The quadratic terms of /j ^ have already been neglected in our pre- 
vious papers. Indeed, studies at the Tevatron suggest that those contribu- 
tions are tiny: See 1 16] for the latest data. 



angular distribution: We can thereby probe exclusively the 
dv term, since any contribution from cIr disappears within 
our approximation as a result of the decoupling theorem 
found in [iMl- 

We are now ready to calculate the following matrices: 



gi(cosef)gj(cos0f) 



Ml, = d cos 9c- — (/, 7=1, 2), 

gi(cos0e) 



m: 



dE, 



hi{Et)hj(Et) 
lii(E() 



(14) 

(/, 7=1,2, 3) (15) 



and their inverse matrices X^j , all of which are apparently 
symmetric. Then the statistical uncertainties for the mea- 
surements of couplings dv and du could be estimated by 



\6dv\ = ^X^^cTiINt = ^Xyi (16) 
through the angular distribution, and 

\6dv\ = ^X^2^t/Ne = ^X^JL, (17) 

\6dR\ = ^X^^crelNt = ^X^JL (18) 

via the energy distribution, where at, Nt and L denote the 
total cross section, the number of events and the integrated 
luminosity for the process, respectively. 

3 Numerical results and discussion 

Below we show the elements of M"^'^ computed for yfs — 
7, 8, 10, and 14 TeV, assuming m, = 173 GeV Q. 

(1) The angular distribution 
(1-1) = 7 TeV 

Mfj = 23.102, = -245.412, = 2607.340. 

(1-2) = 8 TeV 

Mfi = 33.234, = -353.598, M^^ = 3762.753. 
(1-3) = 10 TeV 

Mfi = 59.333, Mjj = -632.179, M^^ = 6736.735. 
(1-4) = 14 TeV 

m;\ = 134.052, Mfj = -1428.300, = 15220.286. 

(2) The energy distribution 

(2-1) = 7 TeV 

Mfj = 23.102, Mfj = -245.412, M'^^ = 0.000, 
Mfj = 2607.658, Mfg = -1.974, M^^ = 15.323. 

(2-2) = 8 TeV 

Mfi = 33.234, = -353.598, Mfj = 0.000, 

Mfj = 3763.252, Mfj = -2.880, M^^ = 21.124. 



(2-3) = 10 TeV 



Mfi = 59.333, Mfj = 



-632.179, Mf,= 0.000, 



= 6737.696, Mfg = -5.120, 



= 35.233. 



(2-4) = 14 TeV 

Mf, = 134.052, Mf, = 



15222.443, M: 



-1428.300, Mfj = 0.000, 



-10.953, 



M: 



33 



72.264. 



Here, all M^^ became zero for the same reason as vanish- 
ing dR terms in the angular distribution, i.e., the decou- 
pling theorem il7l - ll9ll . Using the inverse matrices cal- 
culated from the above elements and Eq.(fT6ll - Eq.lfTSTl. 
we can estimate the statistical uncertainties of the relevant 
couplings. 

Before giving the results, however, we should com- 
ment about instabilities of inverse-matrix computations: 
We noticed that the results fluctuate to a certain extent 
(beyond our expectation) depending on to which decimal 
places of M'^'^ we take into account as our input data. 
Therefore we calculate X''^ for the above M''^ and also 
for what are obtained by rounding those M' '^ off' to two 
decimal places. We then use their mean values to get 6dv 
and 6dR with "errors", which are the differences between 
the mean values and the maximum/minimum. 

Let us explain these treatments more specifically by 
taking (1-1) as an example: For those different inputs, we 
get 



= 1-73 for Mfi = 23.102, ■■■ , 

= 2.57 for Mfi = 23.10, ■•• , 

which lead to the mean value 2.15, the maximum 2.57, and 
the minimum 1.73, and resultant errors +0.42, which are 
derived by 2.57 -2.15 and 1.73-2.15. 

We also have to explain how we can take into account 
QCD higher-order corrections: All the numerical compu- 
tations of Mj^:^ in this section were done with the tree-level 
formulas. In order to include QCD corrections there, we 
multiply the tree cross sections by the ^T-factor. This fac- 
tor disappears in the combination X'::^crf and remains only 
in Ni (= La-() when we estimate 5dvR. Therefore the lumi- 
nosity L discussed in the following should be understood 
as an effective one including K (and also the lepton detec- 
tion efficiency ec). 

We now show the whole results in Table [T] and Table |2] 
for the angular and energy distributions, respectively. The 
uncertainties we encountered here are hard to eliminate, 
but still the results will tell us the necessary luminosity 
for reaching the precision which we aim to realize. For 
example, we need at least L ^ 1500 pb"'(= 1.5 fb"') in 
order to achieve \5dv\ - (9(10^^) in case of measuring the 
angular distribution at the LHC whose colliding energy is 
8 TeV. As mentioned, those L should be divided hy K =^ 
1.5 and ei. If we assume e[ - 0.5, the resultant L increases 
slightly, which however hardly affects our conclusion. 



Table 1. Estimated statistical uncertainties of dy from the 
angular distribution of a decay lepton. 



V? [TeV] 


\6dv\ X VI 


7 


2.15 ±0.42 


8 


2.25 + 0.96 


10 


1.12 + 0.12 


14 


0.73 ± 0.02 



Table 2. Estimated statistical uncertainties of dy and dn from 
the energy distribution of a decay lepton. 



V? [TeV] 


\6dv\ X VI 


\6dR\ X VI 


7 


1.86 ±0.28 


0.35 ± 0.02 


8 


1.70 ±0.51 


0.32 ± 0.05 


10 


0.98 ± 0.08 


0.22 ± 0.01 


14 


0.65 ± 0.02 


0.15 



4 Summary 

Using a scenario of Buchmiiller and Wyler ^ and the 
optimal-observable analysis, we studied the statistical sig- 
nificance of possible non-standard top-gluon couplings 
corresponding to dv,A and top-VK coupling corresponding 
to du as model-independently as possible for the current 
and future LHC experiments. 

Since it has been known that those non-standard cou- 
plings are not so large ['5','6']11 we adopted a linear approx- 
imation for performing the optimal-observable analysis. 
Although the c/a terms disappeared under this approxima- 
tion, we got the following knowledge about the determi- 
nation accuracy of dy and du: 

• Measuring the energy distribution could study both dy 
and dR at the same time with a higher precision than the 
case of the angular distribution. 

• Since the angular distribution is affected by only dy con- 
tributions, this one is suitable for exploring the mecha- 
nism of top-pair productions exclusively though the pre- 
cision of dv is slightly lower than that from the energy 
distribution. 

Furthermore, it was pointed out that there were some 
ambiguities in inverse-matrix computations depending on 
where to round the input data (i.e., M':^-elements) off, 
and consequently estimated statistical errors were affected 
thereby to a certain extent. However, our conclusion about 
the necessary luminosity and new-physics scale A do not 
receive that serious influence, because those ambiguities 
do not change the order of the results. 

Finally, if the LHC is going to be steadily upgraded 
and L = 500 fb"' is achieved at V^ = 14 TeV, dy and 
du could be determined with \6dv\ ~ 0(10"-^) and \6dR\ ~ 
0(10 "*) respectively from the energy distribution: Here, 
\6dv\ ~ 0.001 (\6dR\ ~ 0.0001) means that the contribution 
from effective operator O^^^ (^^liw^ suppressed by A > 

^It is expected that results of LHC experiments with ^fs=S TeV will 
give stronger constraints for dvj^ t2 111 . 



7 (20) TeV if the center value of measured dy (dR) is close 
to zero and C^^, (C^L) ~ 1, which we estimated from 
Eqs.© and (fTTT l. 

On the other hand, concerning c/^ terms, which were 
neglected in this work, they induce CP-violating inter- 
actions, therefore those effects would be better off being 
probed via some kind of asymmetric observables L22il . 
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